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Symmetry properties of the evolution equation and the 
state to be controlled are shown to determine the basic fea- 
tures of the linear control of unstable orbits. In particular, 
the selection of control parameters and their minimal num- 
ber are determined by the irreducible representations of the 
symmetry group of the linearization about the orbit to be con- 
trolled. We use the general results to demonstrate the effect 
of symmetry on the control of two sample physical systems: a 
coupled map lattice and a particle in a symmetric potential. 

PACS number(s): 05.45. -Kb 



I. INTRODUCTION 

Despite the recent wave of interest towards control- 
hng chaotic dynamics ||^,^ an interesting and important 
question of controlling systems with symmetries received 
surprisingly little attention in the physics literature. The 
importance of symmetries in controlling, for instance, 
spatiotemporal chaos is evident, since the systems typi- 
cally show rotational and translational symmetries. Al- 
though the presence of symmetries usually significantly 
simplifies the analysis of system dynamics, it also makes 
control schemes more complicated due to the inherent 
degeneracies of the evolution operators. In fact, the pres- 
ence of symmetries, explicit or implicit, makes a number 
of single-control-parameter methods fail , calling for 
multi-parameter control |^-|^. 

In order to see how these restrictions arise, let us con- 
sider a general discrete-time system (the arguments for 
continuous- time systems are very similar), whose evolu- 



tion is described by the map F : R x R 



R 



N 



,t+l 



(1) 



where z is an A^-dimensional state vector and p is an 
M-dimensional parameter vector. Linearizing about the 
steady state solution of z* = F(z*, p*) and denoting x* — 
z* — z* and u* = p* — p* , we readily obtain 



Ax* +Bu*, 



(2) 



where Aij = dFi (z* , p* ) / dzj is the Jacobian of the trans- 
formation and Bij = dFi{z* ,p*)/dpj is the control ma- 
trix. 

If the steady state z* is unstable, it can be stabilized 
by an appropriate feedback through the time-dependent 
control perturbation u* if the matrices A and B satisfy 



certain conditions. We will understand the design of the 
control scheme as an appropriate choice of the control 
matrix B. We will see below that the conditions affecting 
this choice can be easily obtained from the symmetry 
properties of the system and the controlled state. 



II. STABILIZABLE VS. CONTROLLABLE 
SYSTEMS 

Assuming that the feedback is linear in the deviation 
from the steady state z* we can write 



-Kx* 



(3) 



and one obtains the linearized evolution equation in the 
following form: 



(A-ba:)x*. 



(4) 



The matrix A' — A—BK has stability properties different 
from the stability properties of the matrix A. This can 
be exploited to make the steady state z* (the matrix 
A — BK) stable under control. 

The dynamical system (|^) or the pair (A, B) is said to 
be stabilizable, if there exists a state feedback (||), such 
that the system (^ is stable. The stabiliz ability is the 
property, which often sensitively depends on the values 
of the control parameters p*. 

In the majority of practical situations it is preferable 
to have an adaptive control that would stabilize a given 
steady state z*(p*) for arbitrary values of the system pa- 
rameters. This is especially important if one is to track 
the trajectory z* as p* changes or use the same control 
arrangement to stabilize different steady (or even peri- 
odic) states. 

Such a control scheme is obtained if the more restric- 
tive condition of controllability is imposed on the ma- 
trices A and B. The dynamical system (^) or the pair 
{A, B) is said to be controllable if the eigenvalues of the 
matrix A — BK can be freely assigned (with complex 
ones in conjugate pairs), which is equivalent (see Theo- 
rem 5.13 in Reference [pi) to requiring that: 



rank(C) = N, 



(5) 



where C = {B AB A^B ■ ■ ■ A^'^B) is called the coii- 
troUability matrix. Relation (^ was introduced into the 
physics literature from linear systems theory by Romeiras 
et al. [| as a simple, but practical test of the controlla- 
bility. 
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In order to better understand the restrictions imposed 
by the symmetry, it is beneficial to look at the controlla- 
bility condition written in the form (^) from the geomet- 
rical point of view assuming AI — 1 and B = h. Suppose 
we let the system evolve under control for t steps from 
the initial state x*. The final state will be given by 



m=0 



(6) 



The controllability in this context is equivalent to the 
vectors 



0, 



1 



(7) 



spanning the state space for r = A^, so that any initial 
state can be mapped to any final state in r time steps by 
an appropriate choice of the "coordinates" u*"*""* of the 
difference x*+'^ — ^"^x*. 

If the matrix A is non-degenerate (has a non- 
degenerate spectrum), one can always find a vector b, 
such that the resulting set (0) forms a basis. However, 
if A is degenerate (which is a usual consequence of sym- 
metry), there will exist an invariant subspace L'' C , 
with the dimension dim(L'') > 1, where the dynamics 
of the system cannot be controlled with just one control 
parameter (see Reference for an example of such a 
situation). 

If the system dynamics in happens to be stable, 
the system can still be stabilized similarly to the non- 
degenerate case, but we have to ensure controllability in 
case the dynamics in this subspace is unstable. This can 
be achieved by increasing the number of control param- 
eters M, which extends the set (R), until it spans every 
unstable invariant subspace of R . This would lead one 
to assume that M should be defined by the dimension of 
the largest invariant subspace, or equivalently the highest 
degeneracy of the Jacobian matrix A. We will see how- 
ever, that various kinds of degeneracy have a somewhat 
different effect on the controllability of the system. 



III. THE NUMBER AND SELECTION OF 
PARAMETERS 

Let us assume that the evolution equation (0) possesses 
a symmetry described by a symmetry group ^, i.e. the 
map F commutes with all group actions: 



F{g{z),p) = g{F{z,p)), VgeG, 



(8) 



or in other words, the function F{z, p) is C^-equivariant. 

The symmetry of the linearized equation (||) in the 
absence of control (u = 0) is in general different from 
(although closely related to) the symmetry of the full 
nonlinear equation (^. We will call the respective sym- 
metry group G* . It generates the matrix representation 
T in the state space R^: 



(.g(x)), = (r(g)x), = T{g),,x, 



(9) 



Due to the symmetry all matrices T(g) commute with 
the Jacobian 



T{g)A = AT{g), ^geQ*. 



(10) 



Q* depends on both Q and the reference state z*, or to 
be exact, its symmetry group, which we denote Ti.: 



h{z*) = z*, yhen. 



(11) 



The symmetry of the evolution equation is reduced upon 
linearization, if the reference state has low symmetry, 
and then Q* becomes one of the subgroups of Q. On the 
other hand, Q* might be equal to G or even include Q as 
a subgroup for highly symmetric reference states, with 
the apparent symmetry increased by linearization. 

Decomposing T into a sum of irreducible representa- 
tions of the group Q* with respective dimensionalities 
TTir, r — 1, • ■ • , ^ we obtain: 



(12) 



According to the standard group theoretic analysis g| 
the Jacobian A will have eigenvalues with multiplic- 
ities rir > rrir, corresponding to the dimensions of the 
irreducible representations contained in the decompo- 
sition of T. If rur > 1 for some r, the Jacobian becomes 
degenerate, which causes certain control methods to fail 
(see, for example. Reference [^). 

Next we use the result of linear systems theory based 
on the Jordan decomposition of the Jacobian matrix: 



A = SAS-^ 



A2 



V 



A7 



(13) 



where the Jordan superblock 
/A'-i 

A'' = 

V 



A'- 



A'- 



(14) 



corresponding to the eigenvalue Ar has dimension rir and 
consists of Sr Jordan blocks 



A" = 



1 A. 



V 



A, 
1 



(15) 



Since the controllability is invariant with respect to 
coordinate transformations (Theorem 5.17 in Reference 
H ) , condition (||) is satisfied for the pair (A, B) if and 
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only if it is satisfied for the pair {A,B), where B = SB 
is the transformed control matrix. 

If B is partitioned according to the block structure of 

A: 



B = 



fB^\ 

B^ 






, B" ^ 


\&) 









(16) 



the controllability condition for the pair (A, B) is reduced 
to the controllability conditions for each pair {A^,B'^), 
which, in turn, is satisfied (Theorem 5.18 in Reference 
) if and only if for each r the set of Sr M -dimensional 
row vectors 



irs 



(17) 



is linearly independent. This in turn can be achieved if 
and only if M > for every r. Hence the minimal num- 
ber of control parameters should equal the maximal num- 
ber of Jordan blocks contained in any one superblock: 

= maxsr- (18) 

r 

In general, the system under consideration is not 
Hamiltonian and therefore its Jacobian matrix A is not 
hermitian, hence non-diagonalizable. Therefore we have 
rir > Sr > rur. Howcver in the absence of accidental de- 
generacies Ur = Sr = nir for cvcry r, and the condition 
( p^ ) is equivalent to 



rank(i?r ) 



(19) 



The calculation of the transformation S can be avoided 
for compact groups Q* by using the projection operator 
on the subspace L"" C , that transforms according 
to the r-th irreducible representation: 



P'' = rur 



(20) 



Here x^id) is the character of the group element g in the 
representation and diJ,{g) is the group measure 
For finite groups this integral is replaced with the sum. 
Using the fact that 



rank(B'') rank(B'' 



(21) 



where B^ — P^B, we conclude (assuming there are no ac- 
cidental degeneracies), that the controllability condition 
is satisfied whenever 



M > maxTTir 



(22) 



Therefore the minimal number of independent control 
parameters Mmin is equal to the dimensionality of 
the largest irreducible representation present in the 
decomposition of the representation T of the group Q* in 
the state space . 

Regarding the control matrix _B as a row of M vectors 



B = (bi b2 • • ■ bM), 



(23) 



we see that the control scheme yielding a controllable 
system is obtained by choosing the vectors b^, such that 
TTir of the projections P^hi would be linearly independent 
for every r. 

If accidental degeneracies are present, symmetry prop- 
erties only give a lower bound on the number of required 
control parameters and one should look at the Jordan 
block structure of the Jacobian to determine the control- 
lability using the more general conditions (|l8|) and (p^). 

It is easy to see intuitively why the number of control 
parameters is determined by the number of the Jordan 
blocks Sr and not the multiplicity Ur, if we look at the 
Jacobian already reduced to the Jordan form. For in- 
stance, 



■A 



(24) 



A, 




generates the set of three linearly dependent vectors 
fO b,fi = Ab,f2 = A^b (compare to (0)), that span 
the 1-dimensional subspace of R"^ for any choice of b. 
As a result, 3 linearly independent vectors bi, b2, ba are 
necessary to control the system. 
On the contrary, the Jacobian 



(25) 



generates the linearly independent set of basis vectors 
that spans R^, requiring just one control vector b. 

Finally we should note, that symmetry does not always 
make the Jacobian degenerate, and the non-degenerate 
case can be handled in the same way as the one with 
no symmetries. Neither does the degeneracy by itself 
imply, that multi-parameter control is required. Even 
if > for some r = r' (there is a degeneracy), but 
Sr = mr = 1 for every r (the degeneracy is accidental 
and the Jordan block A^' is not block-diagonalizable) , 
one control parameter is indeed sufficient to make the 
system controllable. 



IV. CONTINUOUS-TIME SYSTEMS AND 
PERIODIC ORBITS 



and nir of M columns of B^ are linearly independent, i.e. 
form a basis in the eigenspace L'' (columns of B^ span 
L'') for every r. 



The obtained results hold for continuous-time systems 
and can be easily generalized to periodic trajectories. We 
should first observe that a periodic trajectory of period 
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T can be treated as a fixed point solution of the super- 
position of T maps. The equation z* = F'^ (z* , p* ) has 
T solutions corresponding to the points of the periodic 
orbit zl, fc = 1, • • • , T. 

Next, we define the time-dependent single-step Jaco- 
bian 



A' 



dzj 



and the control matrix 



dp. 



(26) 



(27) 



The controllability condition for the periodic orbit can 
be generalized requiring that the pairs (Ak,Bk) be con- 
trollable for every k. 

Now, suppose that the symmetry of the state zl is 
described by the group Hk, such that Hk '!= G- We can 
then write 

hizl^,) - h{Fizl)) = FihizD) F(z*) = z*+i (28) 
for every h G Hk, or consequently 



Hi C 7^2 C ... C C Wi, 



(29) 



i.e. the symmetry group of all the states of the trajec- 
tory is the same and can be determined using z^. with an 
arbitrary k: Ti. = Hk- 

This, in turn, means that Q* too is unique for any 
given periodic trajectory as is the representation T. It 
is therefore enough to know the symmetry properties of 
an arbitrary point of the periodic trajectory in order to 
establish the requirements on the control scheme identi- 
cally to the time-independent case. 

And finally, consider a continuous-time evolution equa- 
tion. It can generally be written as 



dtz(t) ^ F{z{t),p) 



(30) 



Linearization of (jS^) around the steady state z* , simi- 
larly to the discrete-time case, yields 



5tx(t) = Ax(t) + Bu{t) 



(31) 



and we again denote Q* as the group of all actions com- 
muting with the action of the Jacobian. The control- 
lability of the pair {A, B) ensures that all eigenvalues 
oi A — BK can be chosen negative, so that the steady 
state becomes stable. As a result, the control matrix B 
should satisfy the same conditions as those obtained for 
the discrete-time case. 



V. COUPLED MAP LATTICE 

Next we apply the general results to the case of the 
coupled map lattice defined by the following evolution 
equation 



= e/(z*_i) + (1 - 2e)/(z*) + e/(4+i), (32) 

with i = 1,2, ■ ■ ■ , N and the periodic boundary condi- 
tions, i.e. 2*_|_jv = imposed. The local map function 
f{x) can be chosen arbitrarily. 

The symmetry group G of the lattice includes trans- 
lations by an integer number of lattice sites (periodic 
boundary conditions make the group finite) and reflec- 
tions about any site. The corresponding point group is 
Cnv It has five irreducible representations T^, - ■ ■ ,T^ 
with mi ~ TO2 = fTT-s = fTT-i = 1 and 777,5 — 2. 

Linearizing eq. (|3|) about the steady state z* we ob- 
tain eq. (|^) with A — CD, where 

a, = (1 - 2e),5,j + e((5,,,_i + 6,,,+i) (33) 

(with Sij±i extended to comply with periodic b.c.) and 

D,,^f{z:)S,,,. (34) 

This partition of the Jacobian into the product of two 
matrices explicitly shows how the symmetry group G* 
depends on the symmetries of the nonlinear evolution 
equation and the controlled state z*. The matrix C has 
all the symmetries imposed by the chosen inter-site cou- 
plings of the nonlinear model: 



T{g)C = CT{g), MgeG, 



(35) 



while the matrix D reflects the symmetries of the refer- 
ence state z*: 



T{h)D = DT{h), yh e H, 



(36) 



and since the Jacobian A only commutes with matrices 
that commute with both C and D, G* should be a max- 
imal subgroup of G and Ti. 

Deriving the restrictions on the control matrix B is the 
next step. Once the group G* is determined, we construct 
its A'^-dimensional representation T and decompose it 
into the sum of the irreducible representations of Cmv 
For instance, a zig-zag state gives G* = Cnv, n — N/2 
and M = — 2; & space-period-s not reflection- 
invariant state corresponds to G* — C'„, 77 = N/ s and 
M — mi = 1, etc. 

In particular, a uniform reference state has G* = Cmv 
and 



(37) 



where each of the representations corresponds to the 
subspace L'^, generated by the Fourier modes (e'^); = 
exp(±zfc^) with wavevectors < |fc| < tt, to fc = 
and T** to fc = tt. Since is present M — 7775 — 2. 
Therefore in order to control an unstable uniform steady 
state of the coupled map lattice we need at least two 
control parameters . This is the reflection of the parity 
symmetry of the model (^2]). 

Choosing B = (bi b2) as a 2-column matrix, and defin- 
ing the Fourier coefficients 
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b'i = e'= • b„ (38) 

we write the conditions on the vectors bi and bi ^ 0, 
^ 0, 6f ^ for < |fc| < TT and either b'l ^ or 
62 7^ for /c = 0, TT. For example, the choice 

Bij = Sj^iSij + (5j,2(5'i,i+i (39) 

yields a controllable system for any 1 < I < N. It corre- 
sponds to applying feedback locally through the pertur- 
bations of the variables at the adjacent sites I and I + 1. 
In fact it can be easily seen that this control arrangement 
makes unstable periodic orbit with arbitrary symmetry 
controllable. 

All the examples above show that the symmetry is re- 
duced upon the linearization of the evolution equation. 
However the symmetry can increase as well. It is quite 
easy to construct a coupled map lattice system whose 
symmetry will increase for certain highly symmetric ref- 
erence states. We will see another (continuous-time) ex- 
ample just below. 

VI. PARTICLE IN A SYMMETRIC POTENTIAL 

The motion of a particle in a symmetric potential, 

mdfr = -yV{r), (40) 

serves as another example of the relation between the 
groups Q and G* ■ Suppose that the potential V{r) pos- 
sesses the cubic symmetry (group O C 50(3)), but is not 
spherically symmetric, for instance: 

V{r) — Vb cosh(fca;) cosh(fc?/) cosh(fcz). (41) 

Linearizing the evolution equation (^ ) about the equi- 
librium point r* = we obtain 

where = —Vok^/m and 1 is a 3 x 3 unit matrix. If 
Vo < the equilibrium is unstable. 

Equation ( |42| ) is spherically symmetric, Q* = 50(3), 
and therefore 5 C tj*, i.e. the symmetry of the linearized 
equation is higher than the symmetry of the original non- 
linear evolution equation. 

Next we notice, that the representation T of Q* in the 
6-dimensional space {r, v} can be decomposed into a sum 
of two 3-dimensional irreducible representations of 50(3) 
(vector representations): 

T^T^®T\ mi =3. (43) 

This indicates that in order to control the unstable state 
r* = V* = one needs at least three independent control 
parameters. 

Probably the simplest way to control such a system is 
re-adjust the potential (applying external fields, shifting 



support point, etc) based on the instantaneous values 
of the position r and velocity v of the particle. This 
corresponds to picking the control matrix in the following 
form: 




where bi,b2,b3 could be chosen as any three linearly 
independent vectors in R^. 



VII. CONCLUSIONS 

Summing up, we conclude that the symmetry prop- 
erties of the system should be understood prior to con- 
structing a control scheme. 

The number of control parameters required to control 
a given state of the system can often be determined us- 
ing only symmetry considerations, without knowing any- 
thing else about the actual evolution equations. The 
knowledge of the evolution equations (at least in the lin- 
earized form) however allows one to choose the control 
parameters (through the matrix B) systematically, avoid- 
ing trial-and-error search. The general idea can be stated 
briefly: the controllability condition requires the control 
arrangement able to break the symmetry of the evolution 
equation completely. 

Care should be taken if there are accidental degenera- 
cies. The knowledge of the multiplicity of the degenerate 
eigenvalues becomes less useful and typically leads to an 
overestimation of the number of control parameters re- 
quired. This case is more complicated and additional 
information about the structure of the Jacobian might 
be necessary in order to determine the minimal number 
of control parameters and construct the control matrix. 
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